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ABSTRACT 

We have investigated the thermal, electrical, and structural properties of Bi and BiCu films that are being 
developed as X-ray absorbers for transition-edge sensor (TES) microcalorimeter arrays for imaging X-ray 
spectroscopy. Bi could be an ideal material for an X-ray absorber due to its high X-ray stopping power and 
low heat capacity, but it has a low thermal conductivity, which can result in position dependence of the pulses 
in the absorber. In order to improve the thermal conductivity, we added Cu layers in between the Bi layers. 
We measured electrical and thermal conductivities of the films around 0.1 K, the operating temperature of 
the TES calorimeter, to examine the films and to determine the optimal thickness of the Cu layer. From the 
electrical conductivity measurements, we found that the Cu is more resistive on the Bi than on a Si substrate. 
Together with an SEM picture of the Bi surface, we concluded that the rough surface of the Bi film makes the 
Cu layer resistive when the Cu layer is not thick enough to fill in the roughness. From the thermal conductivity 
measurements, we determined the thermal diffusion constant to be 2 x 10 3 pm 2 ps _1 in a film that consists of 
2.25 pm of Bi and 0.1 pm of Cu. We measured the position dependence in the film and found that its thermal 
diffusion constant is too low to get good energy resolution, because of the resistive Cu layer and/or possibly a 
very high heat capacity of our Bi films. We show plans to improve the thermal diffusion constant in our BiCu 
absorbers. 
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1. INTRODUCTION 

We are developing a microcalorimeter array for imaging X-ray spectroscopy for future X-ray astrophysics 
missions (Ref. 1). The array has 8x8 pixels. The size of each pixel is 250pm x 250 pm (Figure la). Each pixel 
consists of a Mo and Au bilayer transition-edge sensor (TES) as a thermometer, a silicon-nitride membrane as 
a thermal link, and a Bi and Cu multilayer as an X-ray absorber. We need a close-packed array, so we need the 
absorbers to touch the TES and extend over the space between pixels in a overhanging structure (Figure lb 
and c). 

We make our X-ray absorbers out of Bi because Bi, a semi-metal, has the following advantages over metals. 
Bi has a high X-ray stopping power. Bi has a low heat capacity: in general we need the low heat capacity in 
order to get good energy resolution. Bi does not change the electrical resistance or superconducting transition 
of the TES even when the Bi is directly deposited on the TES. In addition, Bi is brittle, so that it is easy to 
fabricate the overhanging structure with Bi in a micro-machining liftoff process. 
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Figure 1 . Schematic of the absorber, (a) Top view, (b) Side view of old design, (c) Side view of new design. 
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Figure 2. Pictures of the devices used for (a) the electrical conductivity measurements and (b) for the thermal con- 
ductivity measurements. 


However, Bi has a low thermal conductivity, which might result in position dependence of the signals from 
X-rays absorbed within the 250 /max 250 /jm absorber. This would make the energy resolution bad. The 
optimal filter process, which is a standard process to determine the pulse height of the X-ray microcalorimeter, 
assumes that all pulses have the same shape. If the pulses have different shapes, the optimal filter gives us 
wrong pulse heights. Therefore, we put Cu on the Bi in order to make the thermal conductivity higher. Since 
Cu has a high heat capacity, we need to find the optimal thickness of the Cu layer. We would like to match 
the heat capacity to the width of the superconducting transition, so that the dynamic range of the calorimeter 
covers the whole x-ray band. For our typical transitions, we need to augment the Bi heat capacity anyway. 
Ideally, the amount of Cu needed to tune the dynamic range would also provide enough thermal conductivity. 

Considering only the X-ray stopping power and the heat capacity, we have been fabricating absorbers with 
6.5 (im of Bi and 0. 3-1.0 y , m of Cu in total. For most absorbers, we have divided the Bi and Cu into three 
layers each, as in Figure lb. (This multi-layer design was originally implemented to balance stress in thick Bi 
films.) So far, we cannot reject the possibility that bad thermalization in the absorber degrades the energy 
resolution. Since the Cu thickness was chosen based on heat capacity alone, we decided to investigate the 
electrical and thermal conductivities of our BiCu films, using the devices shown in Figure 2. 


2. THEORIES 


2.1. Diffusion Equation 

In order to study the non-steady heat conduction, we need to solve the diffusion equation, 


where T is the temperature, t is time, 
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x, y and z are position, and 
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is the diffusion constant. Here, C is the heat capacity per unit volume and K is the thermal conductivity. The 
propagation time is independent of the cross section but dependent on the square of the distance. 
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Figure 3. SEM pictures of the Bi strips, (a) 2.25 pm thick, (b) 0.25 pm thick. 


2.2. Thermal Conductivity 

If we study the steady heat conduction, although that is beyond the scope of this paper, we can measure the 
thermal conductance of the device, G , using the Fourier’s law, 


G = 


Q_ 

AT’ 


(3) 


where Q is the heat current through unit area per unit time and AT is the difference of the temperature of the 
two thermometers. 


2.3. Wiedemann-Franz Law 

We can estimate the thermal conductivity of the strip from the electrical resistivity, a, which we will measure 
in § 4, using the Wiedemann-Franz law, 

K = LTcr , (4) 

where L = 2.45 x 10 -8 W fi K~ 2 is the Lorentz number. 


3. SEM PICTURES 

We deposited 2.25 pm of Bi on a Si substrate with a deposition rate of 2 nm s -1 by a thermal evaporation 
process without substrate heating. This is our normal deposition method for the Bi, except that we put the 
Bi on a TES that is on a SiN membrane. Figure 3a shows a SEM picture of the 2.25 pm-thick Bi film. The 
Bi film has rough surface with a grain size of ~ 0.1 pm, which is in the same order of the thickness of the Cu 
layer. We also deposited 0.25 pm of Bi with a deposition rate of 0.2 nm s _1 for comparison. We found that the 
thinner Bi film has similar grain size, but has smoother surface than the 2.25 pm-thick Bi film has (Figure 3b). 

4. ELECTRICAL CONDUCTIVITY MEASUREMENTS 

4.1. Setup 

Figure 2a shows the device used for the electrical conductivity measurements. On a Si substrate, we put Au 
bonding pads and deposited a Bi strip on the bonding pads and the Si substrate. On some of the Bi strips, we 
deposited Cu. The length and width of the scrip are 4300 pm and 100 pm, respectively. 

We measured the resistance of the strips as a function of the temperature down to 4 K by 4- wire measure- 
ments, because such measurements enable a quick and comprehensive survey of the BiCu films. 
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Figure 4. (a) Resistance vs. temperature of the Bi strip, BiCu strips and Cu strip, (b) Reproduced resistance vs. 
temperature curves of the BiCu films as a parallel sum of those of the pure Bi film and the pure Cu film. Filled and 
Open symbols represent the actual data and reproduced data, respectively. 


4.2. Resistance vs. Temperature 

Figure 4a shows the resistance vs temperature of a pure Bi strip, four BiCu strips, and a pure Cu strip. 

In both of the thin (0.25 pm) and the thick (2.25 pm) Bi strips, the resistivity increases at low temperatures, 
and their resistivities are similar to each other. In the literature, Bi films also show similar temperature 
dependence of the resistance, although our Bi film has relatively high resistivity in comparison with other Bi 
films (e.g. Ref.2, Ref.3 and Ref.4) suggesting a shorter mean-free path of the electron. 

On the 2.25 /xm-thick Bi strips, we deposited 0.1 ^m, 0.15 pm, 0.2 pm and 0.3 pm of Cu, respectively. At 
low temperatures, since the Bi has high resistivity, the conductivity in the whole strip is dominated by the 
Cu, even in the strip with the thinest Cu layer (0.1 pm ) , while at higher temperatures, the Bi layer makes 
some contribution to the conductivity of the whole strip when the Cu layer is thin (0.1 pm and 0.15 pm). 
Consequently, the BiCu strip with the thin Cu layer shows the temperature dependence similar to the pure Bi 
strips. At 50 mK, we measured the resistance of a BiCu strip with 0.2 pm of Cu and found that the strip has 
the same resistances at 50 mK and 4 K. This result allows us to assume that other BiCu films also have the 
same resistances from 50 mK through 4 K. We use this assumption when we estimate the thermal conductivity 
from the electrical conductivity in § 4.5. 

We measured the resistance of a 0.1 /xm-thick pure Cu strip, for comparison with Cu on Bi. Interestingly, 
the strip shows three times lower resistance than that of the BiCu strip with 0.1 ^xmof Cu. 

4.3. Effective thickness of the Cu layer 

In order to explain the discrepancy in the conductivity for the same Cu thickness, seen in § 4.2, we tried to 
reproduce the resistance vs temperature curves of the BiCu strips assuming that the resistance of the BiCu 
strip is a parallel sum of the resistances of the Bi layer and the Cu layer. We found that we can reproduce the 
BiCu strip with 0.1 pm, 0.15 pm, 0.2 pm and 0.3 pm of Cu as a parallel sum of the 2.25 pm of Bi and 0.03 
pm, 0.06 pm, 0.16 ^m and 0.3 pm of Cu, respectively, as shown in Figure 4b. The results indicate that the 
Cu layer is more resistive on a Bi layer than on a Si substrate, as if the effective thickness of the Cu layer were 
thinner on the Bi layer, and that the effect is more significant when the Cu layer is thinner. 

From these results, together with the SEM picture, we can conclude that the rough surface of the Bi film 
makes the Cu layer more resistive than expected, and that if the Cu layer were thick enough to fill in the rough 
Bi surface and extend as a continuous layer above the Bi, the strip would have the conductivity expected from 
the pure Cu result. The SEM pictures suggest that the effective thickness of a Cu layer could be thicker on a 
thinner Bi film. In the future, we will fabricate such strips to test this expectation. Considering these results, 
we designed a new absorber as shown in Figure lc. We put 1.0 pm of Bi on the TES as an insulator, and 0.5 
pm of Cu on the Bi as a heat bus. We expect that the thin Bi layer and the thick Cu layer make the effective 



thickness of the Cu layer similar to the deposited thickness. On the Cu layer, we put 5.5 pm of Bi as the X-ray 
absorber. 

4.4. Annealing 

It is known that we can improve the electrical conductivity of the Bi and Cu by annealing (e.g. Ref. 4). 
Therefore, we annealed the strips 2 hours at 150 C in vacuum. 

The resistance vs temperature of the Bi strip didn’t change after the annealing. We attempted annealing at 
higher temperatures on a hot plate. However, the Au bonding pads changed color with 15 minutes annealing 
at 200 C and melted with 30 minutes annealing at 230 C, so that we cannot measure the resistance of the 
annealed strips with the Au bonding pads. We are fabricating Bi strips with Cu bonding pads to try the 
annealing at higher temperatures. Note that, however even if we can improve the conductivity, we cannot 
apply the annealing to the absorber that is deposited directly on the Au of the TES. We need to modify the 
fabrication process in order to anneal the absorber on the TES at 200 C or higher. 

The electrical conductivity of the pure Cu strip improved, but the BiCu strips became more resistive due 
to the annealing. Possible explanations of this are inter- diffusion of the Cu into the Bi or Bi Cu alloying, 
although the temperature for the Bi Cu alloying is higher than 150 C. The BiAu strips changed color from gold 
to silver by the annealing, also indicating inter-diffusion or alloying. The rough surface of the Bi supports an 
inter-diffusion hypothesis. 

4.5. Estimation of the Thermal Conductances 

We estimated the thermal conductivity of the Bi strip and Cu strip from the electrical conductivity using the 
Wiedemann-Franz law: K of the Bi is 2.0 x 10 _1 ° W K" 1 ^m _1 and K of the Cu is 1.6 x 10~ 7 W K -1 
K of the Bi is three order of the magnitude higher than that of the Bi. It is known from the literature that C 
of the Bi is 1.5 x 10~ 7 J K _1 cm -3 and C of the Cu is 1.2 x 10~ 5 J K _i cm -3 : C of the Cu is two orders of 
magnitude higher than that of the Bi. Using the C values, we can estimate the diffusion constant of the strips 
as D of the Bi is 1.3 x 10 3 fim 2 /rs _1 and D of the Cu is 1.4 x 10 4 pm 2 ps -1 . Cu has D that is one order of 
magnitude higher than that of the Bi. 

In the horizontal directions, the absorber consists of the Bi layers and Cu layers in parallel, so that we use 
the average K and C of the Bi and Cu in the diffusion equation. If the Cu layer is thick enough to dominate C 
and K of the whole device, D of the whole device approaches D of pure Cu. In other words, there is an upper 
limit on D\ if we use too much Cu, it does not further improve the thermalization. 

In the vertical direction, we have 6.5 /jm of Bi in total without any parallel layer of the Cu. So far, we have 
assumed that we can ignore thermal transport delay in the vertical direction, but we will assess this assumption 
with measurements on pure Bi in the future. 

For the BiCu films, we converted the thermal conductivity into the thermal conductance as 

G = i K = LT b < 5 > 

where A and d are the cross-section and length of the strip, respectively. Here, we used T = 125 mK, and d = 
2180 fx m, which is the operating temperature and the length of the device used for the thermal conductivity 
measurement in § 5, because we will compare the estimated G to the result of the thermal conductivity 
measurement. Filled circles in Figure 4b show the estimated G as a function of the Cu thickness. We also show 
the result from the pure Bi layer in the figure, for comparison. 

So far, we implicitly assumed that the thermal conductance is proportional to the electrical conductance. 
However we note that it is difficult to apply the Wiedemann-Franz law at such a low temperatures. It is known 
that L is much smaller than the value shown above at low temperatures (Ref. 5), so that the estimated thermal 
conductivities are over-estimated. We need the thermal conductivity measurements in order to confirm the 
assumptions. 
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Figure 5. Filled circles show the thermal conductance of the Bi strip and the BiCu strips at 125 mK, estimated from 
the electrical resistivity. A filled star shows measured thermal conductance of the BiCu strip with 0.1 of Cu. 

5. THERMAL CONDUCTIVITY MEASUREMENT 

5.1. Setup 

5.1.1. Devices 

Figure 2b shows the device used for the thermal conductivity measurements. Each device has two TESs and 
a long absorber on a SiN membrane. The membrane has slits along the absorber, so that all the heat in the 
absorber should flow through the two TESs. This device is a kind of Position Sensitive TES microcalorimeter 
(PoST; see Sadleir et al. and Figueroa-Feliciano et al. in this Volume) with a continuous absorber, so that we 
can determine the X-ray hit position without any collimator. The ratio of the pulse heights of the two TESs 
tells us the position. 

The membrane has a thickness of 0.5 /rm. Each TES has a transition temperature of ~127 mK. On the 
membrane and the TESs, we deposited 2.25 /rm of Bi and 0.1 /rm of Cu. The length of the absorber is 2180 
fim, of which 190 /rm overlaps each of the two TESs. The width of the absorber is 100 fim. 

5.1.2. Measurements 

We put the device into an ADR and cooled down the heat bath to 120 mK. We applied bias voltages to the 
two TESs in the PoST. We tuned each bias voltages to obtain similar pulses from the TESs. We irradiated 
the device with characteristic X-rays from Ti. We read out signals from the two TESs using two two-stage 
SQUIDs with a sampling rate of 1 MHz and a cutoff frequency of 300 kHz. We utilized a digital oscilloscope 
to collect the raw data. In order to investigate performance of the devices, we analyzed the raw pulse shapes 
without using an optimal digital filter. 

5.1.3. Event Selection 

We measured the pulse height of each X-ray event as shown in Figure 6a. We see two distributions corresponding 
to the Ti-Ka and /3 lines, although it is not easy to separate them at the ends of the distribution. We also 
measured the area of the pulse for each X-ray event, by integrating the pulse after making the DC level zero. 
We plot the area of the pulse of TES2 against that of TESl in Figure 6b. We see a clear separation between 
Ti-Ka and Ti-K/3 lines. We selected the events in the Ti-Ka line on the area vs area plot for the following 
analysis. 

We divided the X-ray events of the Ti-Ka line into 25 groups, after sorting the events by the area ratio. 
We assumed that each region corresponds to one of the 25 regions of the absorber shown in Figure 2b. We 
calculated the average pulse for each region as shown in Figure 7. 





Figure 6. (a) Distribution of the pulse height of the two TESs. (b) Distribution of the pulse area of the two TESs. 

5.1.4. Monte-Carlo method 

In order to estimate the diffusion constant of the device, we compared the measured result to results of 
simulations. We performed a Monte-Carlo method to solve the diffusion equation, assuming diffusion constant 
of 10 2 , 10 3 and 10 4 /rm 2 /rs -1 . We assumed a 10 p m x 10 jum mesh. We considered only the horizontal diffusion, 
because the thickness of the film is negligible. 

5.2. Diffusion Constant 

5.2.1. Position Dependences 

We measured the pulse height, the pulse area, and the rise time of the 25 average pulses as a function of 
position as shown in Figure 8. We see a saturation of the pulse height in Figure 8a. This is because the pulse 
was not small enough to keep the TES resistance on the transition edge, due to too little heat capacity. We 
need the low heat capacity for the good energy resolution, but if it is too low, we suffer from the saturation 
of the pulses. We also see the saturation of the pulse area in Figure 8b, but the saturation is less significant. 
This is the reason why we used the pulse area ratio instead of the pulse height ratio to determine the position. 

We superposed the Monte-Carlo results on Figure 8 as open circles, open squares, and open triangles. We 
found that diffusion constant of 2 x 10 3 pm 2 ps _1 reproduces the measured data well. 

5.2.2. Delay of rise 

Figure 9a shows average pulses of the TESl and TES2 when X-rays hit the absorber on the TESl. The pulses 
of the TESl and TES2 were collected by the oscilloscope at the same time. We see a delay of the rise in the 
TES2 pulse. This means that it takes ~100 ^is for the heat to travel 2180 /r m. Figure 9b shows the Monte-Carlo 
pulses of the TESl and TES2 for the same X-ray hit position that generated assuming a diffusion constant of 
10 3 pm 2 /rs _1 , whose delay of rise is in good agreement with the actual data. 

5.2.3. Theoretical value 

We estimated the diffusion constant assuming the free-electron approximation, 

D=^v ( l, (6) 

where Vf is the Fermi velocity and l is the mean free path. If we assume the inferred effective thickness of 
the Cu, 0.03 jam, as the mean free path, we obtain 1.6xl0 4 fj, m 2 /is _1 . This is one order of magnitude larger 
than the actual diffusion constant. This means that we cannot apply the free-electron approximation at such 
a low temperature, so we need the actual measurement of the thermal conductivity to estimate the diffusion 
constant. 





Figure 7. Average pulses of TES1. Each pulse corresponds to one of the 25 X-ray hit positions. 


5.2.4. Comparison with the electrical conductivity results 

In Figure 4b, we superposed the actual G determined from the thermal conductivity measurement as a filled 
star. We converted the diffusion constant into the thermal conductivity using Eq. (2). If we consider the overes- 
timation in the Wiedemann-Franz results, the results of the electrical and thermal conductivity measurements 
are consistent with each other. 

5.3. Energy Resolution 

We estimated the degradation of the energy resolution due to the position dependence. 

5.3.1. Energy loss 

In the PoST, the X-ray energy is divided into the two TESs, so that the pulse area has the position dependence. 
However, this does not mean that the pulse area is dependent on the X-ray hit position in the standard TES 
calorimeter. In the area vs area plot, the events are on straight lines. Since the area of the calorimeter pulse 
is proportional to the X-ray energy, this means that there is no energy loss, namely all the X-ray energy is 
collected by the two TESs. In the standard TES calorimeter, there is only one TES, so that all X-ray energy 
is collected by the one TES, namely the pulse area should be independent of the X-ray hit position. 

5.3.2. Pulse shape 

In Figure 10, we show average pulses of two adjacent regions, after correcting the pulse scale to give equal area, 
considering the result in § 5.3.1. The difference of the pulse heights in the two regions is ~20% of the pulse 
heights. This means the thermal diffusion constant of the film is too low to get the good energy resolution. 
Fortunately, 0.1pm of Cu is much thinner than the amount needed to meet the design heat capacity in the 250 
pm wide pixels, for typical transition widths and using published values for heat capacity in Bi and Cu. 

There is a possibility that our Bi film has a much higher heat capacity than expected, so that the 0.1pm 
Cu layer is not thick enough to dominate the heat capacitance in the whole absorber. We can determine the 
thermal conductivity of the SiN membrane by the IV measurements, namely by applying a bias voltage to 
the device and measuring the current through the TES. We can roughly determine the ratio between the heat 
capacity of the device and the thermal conductance of the SiN membrane using the decay time of the pulse. 
Consequently, we can determine the heat capacitance of the device using the IV curves and the decay time of 
the pulse. In the present case, we obtained a heat capacitance that is 1-2 orders of magnitude higher than the 
expected value. We need to understand the Bi contribution to the heat capacity before we can determine the 
optimal Cu thickness. 
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Figure 8. Position dependences of (a) the pulse height, (b) the pulse area and (c) the rise time, (d) Rise time vs pulse 
height. Filled circles represent the actual data. Open circles, open squares and open triangles represent the simulated 
data with a diffusion constant of 10 2 , 10 3 and 10 4 ^m 2 /rs _1 , respectively. 
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Figure 9. Comparison of the TES1 (dotted line) and TES2 (solid line) pulses when an X-ray hits the absorber on the 
TESl. (a): actual pulses, (b): simulated pulses. 
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Figure 10. Comparison of the pulse shape of the 20th (solid line) and 21st (dot line) regions of the 25 regions, shown 
after the area correction. 


6. FUTURE PLANS 

6.1. Improvement of the Bi films 

We will try to improve our Bi film, especially the smoothness of the Bi surface, by annealing after the deposition, 
by heating the substrate during the deposition, and/or by slow-rate deposition. We will test these devices by 
taking the SEM pictures and by measuring the resistance versus temperature. 

6.2. Electrical Conductivity Measurements 

We will also deposit Cu on the thinner Bi and measure its electrical resistance to confirm the expectation that 
the rough surface of the Bi film makes the Cu resistive. 

6.3. Thermal Conductivity Measurements 

We will fabricate and test a device with pure Bi to measure the thermal conductivity of the Bi film itself and 
examine whether the vertical thermal resistance in the Bi is negligible or not. If it is not negligible, we need to 
improve the design of the absorber, for example by adding a vertical Cu layer for the thermal conductivity. 

We will fabricate and test the thermal conductivity measurement devices with thinner Bi and/or thicker 
Cu. If the resistivity of the Cu is improved thanks to the use of a single thicker Cu layer and the smoother 
surface of the thin Bi, the thermal conductivity of the Cu increases but the heat capacity of the Cu stays the 
previous value, so that we can get a larger diffusion constant. We will also fabricate more thicker Cu devices 
to test whether the diffusion constant becomes higher or not. 

In this paper, we investigated the time-dependent heat conduction using the X-ray pulse as a heat source 
and measured the diffusion constant. Added to this, if we can study the steady heat conduction to measure the 
thermal conductivity directly, we can determine C of our Bi film, which is suspected to be much higher than 
the normal value, by dividing the diffusion constant by the thermal conductivity. We can use the two TESs of 
the device as the thermometers for the thermal conductivity measurements and can also use one of them as a 
heater to supply the heat current. 
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